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Abstract. The much wider transverse-momentum range accessible in heavy-ion 
colUsions at RHIC and at the LHC allows us to disentangle the dynamics of partonic 
equilibration from the dynamics of delayed hadronization. This provides a novel tool 
for testing the equilibration mechanisms underlying QCD thermodynamics. Here, I 
argue, on the basis of simple formation-time arguments, why this is so, and I review 
recent theoretical developments in this context. 
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1. Hadronization vs. thermalization 

How can one determine, on the basis of the measured hadronic final state, whether 
equihbration processes occur in a medium of rapidly decreasing density, and if so, 
whether they occur in terms of partonic or hadronic degrees of freedom? The wide 
transverse-momentum range accessible in nucleus-nucleus collisions at collider energies 
provides a novel access to this long-standing question. 

To see this, consider a parton of high transverse energy Et, produced in some hard 
collision: if the parton escapes into the vacuum, then it will reduce its initial virtuality Q 
by perturbative parton splitting. After some time ~ 1/Qhadr, its virtuality is degraded 
to a hadronic scale Qhadr ~ 1 GeV. Hadronization is the non-perturbative dynamics of 
the further fragmentation of this multiparton object. Numerical estimates for the time 
scale of hadronization vary significantly [HEIISI, but owing to the Lorentz boost to the 
laboratory frame, they are proportional to the energy, Lhadr ~ 0(1)-tt^ — tt^- 

What happens if the hard parton escapes into an infinitely extended quark gluon- 
plasma instead? Because of medium-induced gluon radiation, the initial perturbative 
parton splitting is even more efficient. However, the parton cannot hadronize in the 
dense medium. Instead, after some time, its partonic fragments can no longer be 
distinguished from the heat bath: the hard parton is thermalized. To estimate the time 
scale Ltherm for this process, require that the hard parton has lost all its energy through 
medium-induced gluon radiation. According to the BDMPS energy loss formula |3], 
Erp ~ AE = ^x^^-^therm- The partonic thermalization length is Ltherm ~ V^t- 

The above time estimates are simplified and may be improved They illustrate, 
however, that for large transverse energies Ej-, perturbative equilibration mechanisms 
can remain undisturbed by hadronization over a significant time scale, see Fig. ^ 
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Figure 1. Comparison of the hadronization time scale ihadr = n n 

the partonic thermalization time scale Lthcrm = \/^/o:sCf\/ Et/q- For hadronization, 
Qhadr = 1 GeV and c = 2 to account for multiple parton branching. For thermalization, 
otsCp = 1 and g = 1 GeV^/fm. See text for motivation of parameter values from data. 
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Depending on its in-medium pathlength Lmed, the hard parton will either be absorbed 
(-^therm < -^mcd < -^hadr), OT it has a Sufficiently large transverse energy to suffer 
only the onset of equilibration processes {L^ed < -^^therm < -Z^hadr)- In the latter case, 
the parton appears as a medium-modified jet. For lower transverse energies, there is 
not only a competition between the hadronization and the thermalization mechanism 
(-^hadr ~ -^therm)- There is also the possibility that the medium interferes with the 
dynamics of the hadronization process (Lhadr ~ -^med)- For even lower transverse 
momentum, the hadronization time scale is determined by the density evolution of 
the medium, which is not accounted for in the above estimates. Only in this kinematic 
"bulk" regime may formed hadrons stay in contact with the equilibrating medium for 
a significant duration. I now discuss in more detail the recent theoretical developments 
in these different kinematic regimes. 

2. Hard pt 

For sufficiently hard partons, the hadronization time scale exceeds the partonic 
thermalization time scale, Lhadr > -^therm- The current understanding of suppressed 
leading hadroproduction in nucleus-nucleus collisions is based on this inequality. One 
starts from the assumption that the hard parent partons are produced perturbatively, i.e. 
proportional to the number of binary collisions. The production rate of their daughter 
hadrons, however, is reduced, because some parent partons are absorbed in the medium 
(Lined > -^therm) whereas the others {L^ed < -^therm) lose an additional fraction of their 
initial energy through medium effects and thus materialise in hadrons of lower transverse 
momentum. As a consequence of this picture (see Fig. |2I), all observed high-pj- hadrons 
have hadronized outside the medium. Since hadronization of a parton depends only 
on its identity but not on its history. Fig. |21 implies that the ratio of different identified 
leading hadrons stays unchanged with respect to p-p collisions while their absolute yield 
is suppressed with respect to the binary scaling assumption. In Au+Au collisions at 
RHIC, this is observed in the particle species dependence of the nuclear modification 
factor jElIZI, which ceases for hadronic transverse momenta p^'^'^ > 6-7 GeV [Zj (this 
corresponds to higher partonic momenta p^^^°^ > 9-10 GeV). This indicates that for 
the observed hadrons, Lhadr ~ -^therm < -^med at p^^^°'^ ~ 10 GeV. This scale compares 
well with the estimate in Fig. 

To account for the partonic medium-modification prior to hadronization, we must 
understand how the parton showers associated to hard partonic production processes 
interact with a dense medium. This is a difficult problem, which is not solved completely. 
All existing approaches [1101211101111] indicate that the leading parton loses additional 
energy owing to medium-induced gluon radiation and that the parton shower broadens 
in transverse momentum. Remarkably, only one property of the medium enters these 
calculations, namely the amount q of transverse momentum squared that the medium 
transfers to the parton per unit pathlength. Momentum broadening is characterised 
by {k\) ~ qL^ed |4) the additional radiated energy is determined by the characteristic 
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Figure 2. Schematic view of medium- modified leading hadroproduction: the 
hadronization process is assumed to be separated in space and time from the medium- 
modification of the partonic evolution. This picture implies a trivial particle species 
dependence of the nuclear modification factor and applies at RHIC above p^'^^°'^ ~ 10 
GeV only. 



energy scale Uc = and medium effects regulate the additional gluon radiation 

in the infrared on a scale uj ~ uJc/ {(jL^ed)'^^^ [l^J- Also, the expansion of the medium 
can be controlled in these calculations ^3 |^ . Since the transport coefficient is 
proportional to the density of scattering centres, it decreases in time, (?(t). The resulting 
medium-induced radiation spectrum corresponds to the spectrum of an equivalent 
static medium in which all kinematic variables are rescaled by the linear line-averaged 
transport coefficient q = jj— ^^'^"'^ t qir) dr [THj. In practice, this dynamical scaling 



law allows a simplified data analysis in terms of static medium properties; a posteriori, 
one then translates the extracted static average transport coefficient into the realistic 
dynamical one. 

So far, comparisons of this formalism to data are limited to the leading hadron 
spectra and leading back-to-back correlations. In this case, the reduced energy of the 
leading hadron directly traces the reduced energy of the leading parton (see Fig. 12)); 
the reduction is expected to depend mainly on the parameter combination qLi^^^. 
The most advanced comparisons with data [ Q reproduce the absolute scale of the 
nuclear modification factor and its centrality dependence. This allows us to disentangle 
information on geometry (Lmed) and density (g). A satisfactory agreement with data 
is found for an initial energy density of el^-o = 15GeV/fm'^ at initial time tq = 0.2 
fm/c jl]. In contrast to other determinations, this is an energy density with which a 
hard parton was interacting. The value is consistent with the line-averaged transport 
coefficient q ~ IGeV^/fm used for the estimate in Fig. [T] Based on ej^-o, one can 
estimate the duration tqgp over which the density was above the critical energy density 
ec ~ 0.75 GeV/fm^ predicted by lattice QCD. One finds tqgp = 4±2±? fm/c, where the 
central value is for a one-dimensional Bjorken expansion (tqgp = Toero/ec) and ±2 fm/c 
is associated to model-intrinsic uncertainties. These stem from our limited knowledge of 
the dynamical expansion and evolving geometrical distribution of scattering centres, and 
from significant uncertainties related to the high-energy approximations employed by all 
current calculations (for a quantitative discussion, see Refs. Other systematic 
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errors (e.g. those associated to the perturbative calculation of medium-induced gluon 
radiation) are more difficult to quantify. 

How can we further test the partonic mechanism of medium-induced gluon radiation 
assumed to underlie medium-modified high-p^- hadroproduction? Current theory 
predicts i) the dependence of medium effects on Lmed, ii) their dependence on the 
parton identity, and iii) the relation between leading parton energy loss and transverse- 
momentum broadening of the parton shower. In my view, the perspectives for refined 
tests are the following: 

i) Dependence on geometry: The approximately quadratic dependence on the in-medium 
pathlength Lmed is a characteristic feature of non-abelian partonic energy loss jl]. 
Typical observables sensitive to it are the azimuthal anisotropy V2{pt) at high px, the 
centrality dependence of leading hadron spectra and their back-to-back correlations. In 
practice, however, uncertainties in modelling the geometry and dynamics of the collision, 
as well as current experimental systematic errors make it difficult so far to discriminate 
even between a linear and a quadratic dependences on L^aed P^ . 

ii) Dependence on parton identity: a) Hard gluons are expected to suffer larger final- 
state medium effects than hard quarks since their coupling is stronger. Sensitive 
to this effect is, for example, the ratio of any pair of leading hadrons that receive 
different contributions from parent quarks and parent gluons (e.g. the p/p ratio). In 
practice, however, an analysis of these measurements is complicated by the experimental 
limitations on high-pT- particle identification, and by the uncertainties in the vacuum 
fragmentation functions at large momentum fraction z ^Hj- At sufficiently high 
luminosity and/or centre-of-mass energy, a more promising test of gluon propagation 
in matter may be 7-jet or three-jet [qqg) events where tagged leptonic or hadronic 
high-pT- information signals the presence of a gluon. b) Massive quarks are expected to 
show smaller medium modifications than massless ones and this should show up in the 
leading charmed and beauty hadrons (201 El EH 122] • However, for massive quarks, Lhadr 
is reduced since massive quarks are slower, and Lthorm is increased since they lose less 
energy. As a consequence (see Fig. H}, the transverse energy scale where Lhadr ~ -Z^therm 
and partonic energy loss calculations start to apply may lie at the upper end of the 
current experimental data, above p^'^'^ ~ 10 GeV. 

iii) Beyond leading hadron spectra: Calculations predict how the energy lost by the 
leading parton is redistributed in transverse phase space [23 II2|- Ideally, this can 
be tested by reconstructing the total energy associated to the initial parent parton, 
i.e. by measuring the "medium-modified jet". The corresponding jet production cross 
section is expected to follow binary scaling, but the jet shape is broadened and the 
jet multiplicity is softened and increased j21]. The main problem of calorimetric jet 
measurements is the underlying high-multiplicity background. At present, this is studied 
intensively for the heavy-ion programme at the LHC, and first jet observables are known, 
which are sensitive to transverse-momentum broadening but insensitive to the low-prp 
background j21] . Where calorimetric measurements are not feasible, a pragmatic widely 
adopted alternative is the study of jet-like near-side particle correlations associated to 
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high-pT- trigger particles 123- Finding in this way the lost remnants of jets is certainly an 
important qualitative support of the assumption of medium-induced final-state energy 
loss. However, it is difficult to compare the trigger selection of a special subclass of all jets 
(namely those with a particularly energetic leading hadron) with current parton energy 
loss calculations without significant further assumptions. This is so because theory is 
more reliable for energy distributions than for multiplicity distributions (which require 
a more detailed knowledge of hadronization) , and because the same leading-particle 
trigger selects significantly different jet classes with different average Et depending on 
the a priori unknown medium modification. 
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Figure 3. Schematic view of the dynamics of hadronization. In nucleus-nucleus 
collisions, the medium allows a test of the non-perturbative long-range part of parton 
fragmentation if the transverse momentum is not too large. 



3. Intermediate px 

Hadronization is the non-perturbative process that bleaches coloured partons into sets of 
hadronic states. As of today, its dynamics remains a black box (see Fig.|2I) which is filled 
differently by different hadronization models. A satisfactory theoretical understanding 
exists only at sufficiently high pt when factorization theorems are applicable. This black 
box can then be bypassed by a simple mapping, the fragmentation function Dh/q{z, Q"^). 
In this case, the only dynamical information accessible about hadronization is the 
dependence on the scale at which perturbative physics is interfaced with Dh/q{z, Q"^). 
This hmits the possibility to study how hadronization proceeds. 

The same mapping Dh/q{z,Q'^) is applied in the current description of medium- 
modified hard processes (see Fig. EI). As a consequence, the particle ratio of leading 
protons/pions is predicted to trace the ratio of the corresponding fragmentation 
functions, p/n ~ Dp/g{z,Q'^)/DT^/g{z,Q'^) ~ 0.1. However, in Au+Au colhsions for 
phadron ^ ^ QcV, ouc observes a marked deviation [7j from this perturbative ratio, 
which increases with centrality and reaches unity. Thus, perturbative fragmentation 
breaks down in Au+Au collisions although it is known to provide a satisfactory 
description of the particle species dependence of elementary e~^e~ or hadronic collisions 
in the same intermediate-pr regime. The scale p^'^^ < 7 GeV where this happens lies 
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in the range where Lmcd ~ -^hadr ~ -^thorm (see Fig. H)). This suggests that the much 
larger spatial extension of dense matter in nucleus-nucleus collisions interferes with the 
hadronization process at higher pt- The medium becomes a tool for testing the dynamics 
of hadronization. From this perspective, the theoretical challenge at intermediate pt is 
to discriminate, e.g. between the different proposed prehadronic objects through which 
hadronization may occur. 

At present, such a program has at best started. Different models, which focus on 
hadronic spectra at intermediate pt i) interpolate between soft hydrodynamic evolution 
and hard high-pj- processes [2^], including known nuclear effects h) model the 
dynamics of string fragmentation |2ZI, or iii) aim to account for data in terms of the 
recombination of dressed quarks (2H1 1211 IHOl IHI] • In all these approaches, the particular 
non-perturbative model components become numerically unimportant outside the 
intermediate-pT regime above p^'^'^ ~ 7 GeV [2111 • In particular, recombination models 
offer an explanation for the anomalous particle species dependence at intermediate pt 
in terms of a simple prehadronic picture. For example, they emphasise that the elliptic 
flow of different baryons and mesons falls on an universal curve \v2{pT/n) if rescaled by 
the number n of valence quarks pB. • This may be interpreted as a remnant of a common 
underlying partonic flow, which recombines differently for different hadron spectra |32j . 
It opens the exciting possibility that the partonic flow saturates at a rather small value 
of ~ 7% but that it is amplified by a factor 2-3 in the hadronization process. Whether 
this suggestive counting rule can result from a dynamical context (which should explain, 
for instance, the fate of gluonic degrees of freedom) and how it can be discriminated, 
e.g. from models of string fragmentation, remains to be established. 

In short, the strongest argument for the novel physics potential at intermediate px 
is the coincidence of scales Lmed ~ -^therm ~ -^hadr, which indicates the interference of 
the hadronization process with the medium. The data show clear signals of such an 
interference [EllZI- 

4. Soft Pt 

Hadronic spectra at high and intermediate pt differ significantly from thermal and chem- 
ical equilibrium distributions, but theory relates their medium-dependence dynamically 
to the onset of partonic equilibration and medium-dependent hadronization. At soft 
Pt {pt < 2 GeV), the challenge is the inverse: many bulk observables show (approxi- 
mate) equilibrium distributions and the task is to establish to what extent this is the 
consequence of equilibration processes. 

In particular, particle ratios do not show significant deviations from chemical 
equilibrium and can be described by the model of statistical hadronization in terms of a 
temperature (Tch ~ 170 MeV) and a baryochemical potential (/xb ~ 40 MeV at ^/snn = 
200 GeV) jnni- However, the observed equilibrium does not automatically imply a 
dynamical equilibration mechanism, since statistical particle production according to 
the principle of maximum entropy is sufficient to motivate the underlying model. A 
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dynamical origin of the observed equilibrium distribution still remains an intriguing 
possibility in particular because the extracted temperature and baryochemical 
potential lie close to the phase-transition line of equilibrated lattice QCD [HSl EEl EZl ESj • 
But to establish a firmer link between the phase-transition line to a QGP and the 
statistical description of hadron yields is most likely to require a better understanding 
of the hadronization process. 

Moreover, combinations of the model of statistical hadronization with hydrody- 
namically motivated parametrisations of the bulk dynamics at kinetic freeze-out (the 
so-called "blastwave models" ) account well for the particle-species dependence of soft pt 
hadronic spectra and their azimuthal dependence V2{pt) ,39 . [Intriguingly, this implies 
that at soft px the species dependence of elliptic flow is a mass effect, while it is argued 
to be a baryon/meson effect at intermediate pt jSHj-] These combinations also account 
for the space-time structure reflected in identical particle correlations, and in particular 
for the recently measured azimuthal dependence of these correlations, which trace the 
space-time picture of elliptic flow. The model parameters extracted from fits to data 
suggest that the kinetic freeze-out temperature lies significantly below the chemical one 
(Tch — Tkin ~ 70 MeV ~ O(Tkin)) and that the system shows a strong collective trans- 
verse expansion vt > \c. They also suggest a rather short freeze-out time of r/ < 10 
fm/c. [This bound is based on the assumption of Bjorken longitudinal expansion. If the 
longitudinal expansion is slower as, for instance, suggested by the description of dN/ dy 
in terms of Landau hydrodynamics jlO], the freeze-out time can then be significantly 
larger, see Fig. 1 of jH].] However, blastwave models are only (physically motivated) 
parametrisations of data. They do not address the question of which microscopic dy- 
namics leads to the observed distributions. 

Hydrodynamics, based on the assumption of perfect local thermal and chemical 
equilibrium, is a dynamical explanation. It can account for hadronic transverse [I^ H5] 
and longitudinal |12j single-particle spectra and their azimuthal dependence 
In contrast, hydrodynamic calculations do not reproduce the identical two-particle 
correlations jlH]. It remains unclear whether this is a problem of hydrodynamics 
itself |42], or of the freeze-out condition with which hydrodynamics is interfaced [IHl 
Wf\ EHj. The reason is that HBT measurements test the surface of last scattering, 
which can change significantly depending on the freeze-out condition while the large 
mean-free path makes hydrodynamical assumptions questionable. On the other hand, 
a clearly model-intrinsic problem of hydrodynamics is to understand microscopically 
how a hydrodynamic behaviour comes about [IH]- Simulations indicate that for the 
parton densities and partonic cross sections expected at RHIC, the pressure built up 
in a peripheral collision can result in an elliptic flow of ~ 5% at most. This falls short 
of the ~ 15% observed on the hadronic level. Here, the novel provocative suggestion 
from recombination models is that this small partonic V2 may be almost sufficient since 
it gets amplified in the hadronization process 13^. It is important to establish whether 
this is a viable explanation, since this would also strongly affect the understanding 
of V2 ai Pt < 2 GeV. If true, it implies that insight into the hadronization process is 



Theoretical Overview 



9 



indispensable for understanding whether and how hydrodynamics works on the partonic 
level. 

5. Saturation physics - A complementary view 

The discussion of Fig. ^ is not the only way to attach fundamental physics questions to 
the transverse momentum scale. Saturation physics provides an alternative view. Its 
starting point is the observed growth of parton distribution functions at small x. On 
general grounds [SII, this growth must be tamed eventually by density effects and it 
will saturate for densities of order ~ l/a^- In this regime, perturbation theory in a 
high-density background allows for the derivation of density-modifications to the small- 
X BFKL evolution |nilE2lEn]- As a consequence, the unintegrated gluon distribution 
is depleted (saturated) up to a saturation scale Qs jHH (HHl IHE] • This scale can grow 
perturbatively large (possibly Ql ~ (l-2GeV)^ at RHIC), and it increases with the 
nuclear extension (Q^ ~ A^^^) and for small x or increasing rapidity y (Qsiv) ~ Ae"^*^^, 
c ~ 0(1)). In addition to the perturbative high-p^ region in which density effects 
are negligible, saturation physics identifies two transverse-momentum regimes with 
qualitatively different physics: 

Soft pt- saturation region. Since the initial parton density is depleted, one may 
expect that less partons will be produced in the colhsion. It has been argued that this can 
explain the low event multiplicities at RHIC [SZI and that it provides a fair description 
of bulk observables such as rapidity distributions and their centrality dependence jSHl- 
However, a significant mo del- dependent input is required to interface saturation physics 
with soft bulk observables, and models based on other equally fundamental physics 
assumptions are equally successful, see e.g. |I2] . 

Intermediate pt '■ scaling window. The low-fcy (^t < Ql) gluons of a depleted gluon 
distribution are the dominant source of high-Zc^ gluon production of during the small- 
X evolution. Hence, small-x evolution reduces the growth of the gluon distribution 
above the saturation scale, in the so-called scaling window Q^(xevoived) < < 
Qt {x evolved) / Qlix initial) [S^ EUl Elj • Siuce the Saturation scale grows with x, this scaling 
window opens as a function of increasing cms energy or rapidity. One therefore expects 
in both hadron-nucleus and nucleus-nucleus collisions a suppressed hadroproduction at 
intermediate pr iff the saturated gluon distribution is evolved sufficiently far in x jH^ . 
At mid-rapidity at RHIC, this scenario is ruled out by the observed opposite centrality 
dependence of the nuclear modification factor in d-Au and Au-Au [HBl lEH ^Ei IHH] . 
This does not exclude the parton saturation as a viable explanation in the soft regime 
Pt < Ql RHIC mid rapidity. It just shows that the saturated gluon distribution was 
not yet evolved sufficiently far in x for the scaling window to open. 

The observed disappearance of the Cronin enhancement in d-Au at forward 
rapidity [HII was predicted as a consequence of the small-x evolution of a saturated 
gluon distribution [Hllinni- In the context of saturation physics, this signals that the 
scaling window has opened because of the small-x evolution. However, the agreement 
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is at best qualitative: the studies jSHl EH IHS] calculate, on the basis of a kinematically 
simplified partonic cross section, a ratio of gluon spectra; neither the validity of this 
simplification, nor the effect of neglecting quarks has been discussed so far. Also, 
alternative explanations are not fully explored yet. In general, if multiple scattering 
leads to an effective redistribution of hadronic yields not only in transverse but also 
in longitudinal momentum, this could then deplete the ratio -RdAu at forward rapidity 
equally well. Calculations of -RdAu based on multiple scattering [701 El did not confirm 
this effect, but predicted a Cronin enhancement at forward rapidity instead. However, 
these calculations are based on a high-energy approximation in which a transverse- 
momentum transfer does not degrade longitudinal momentum and where only elastic 
multiple-scattering contributions are considered. Thus, it is still conceivable that 
refinements change the result of conventional Cronin calculations significantly. 

If the interpretation of the rapidity dependence of -RdAu in terms of saturated gluon 
distributions stands further tests, then this has far-reaching consequences. At RHIC, 
one expects further observable traces, in particular for open charm production [Z2l IZH] 
and possibly for dileptons 74, . For the LHC, a confirmed onset of parton saturation 
at forward RHIC rapidity implies that the medium attained in Pb-Pb collisions at the 
LHC differs qualitatively from that at RHIC, since it corresponds to largely evolved 
saturated parton densities. In the long term, this finding would give further support to 
the already strong case for a dedicated e-A experimental program. 

6. Conclusion 

In summary, the heavy-ion programs at RHIC and the LHC are in the fortunate situation 
that hadronization and partonic thermalization length scales become comparable in the 
very range of in-medium pathlength and transverse momentum accessible at collider 
energies (see Fig. P). This suggests that the much wider p^-range at RHIC and LHC 
provides a novel tool to disentangle the two recurrent dynamical features of heavy- 
ion collisions: the dynamics of partonic thermalization and the dynamics of medium- 
dependent hadronization. I argued why the main open questions at soft pt are linked to 
our poor understanding of hadronization and freeze-out, and how the study of data 
at intermediate px may improve this understanding. At high px, I discussed the 
possibility to study how an out-of-equilibrium parton - undisturbed by hadronization 
-evolves towards partonic equilibrium. This is not only an important step towards QCD 
thermodynamics. A refined theoretical understanding of medium-modified hard probes 
also allows for more detailed tests of the dense bulk matter with strongly interacting 
penetrating probes. 

Clearly, the RHIC-ness of scales accessible at collider energies is not exhausted by 
this discussion. Varying other scales such as the ^-dependence and energy dependence 
of nucleus-nucleus collisions provides complementary information, for example for a 
precise scan of energy density and baryochemical potential. Also, other views on the 
wide Pt range, such as saturation physics, may turn out to be more fruitful in the 
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future. Whatever this future may be, progress in the strong interactions has always 
been based on the interplay of theory and experiment. The number of accessible scales 
and fundamental concepts illustrates clearly that at present one of the most active and 
most diverse fields for this interplay is the study of nucleus-nucleus collisions. 
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